Interaction Effects in Conductivity of Si Inversion Layers at Intermediate 

Temperatures 
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We compare the temperature dependence of resistivity p(T) of Si MOSFETs with the recent 
theory by Zala et al. This comparison does not involve any fitting parameters: the effective mass 
m* and <?*-factor for mobile electrons have been found independently. An anomalous increase of p 
with temperature, which has been considered a signature of the "metallic" state, can be described 
quantitatively by the interaction effects in the ballistic regime. The in-plane magnetoresistance p(Bu ) 
is qualitatively consistent with the theory; however, the lack of quantitative agreement indicates that 
the magnetoresistance is more susceptible to the sample-specific effects than p(T). 
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The theory of quantum corrections due to single- 
particle (weak localization) and interaction effects has 
been very successful in describing the low-temperature 
electron transport in low-dimensional conductors (see, 
e.g., [I]]). There was, however, a noticeable exception: 
the temperature- and magnetic-field dependences of the 
resistivity p of the most ubiquitous two-dimensional (2D) 
system, electrons in silicon MOSFETs, defied the the- 
oretical predictions. Serious quantitative discrepancies 
between the experimental data and the theory of inter- 
action corrections were noticed two decades ago (see, e.g., 
The gap between the expected and observed low- 
temperature behavior of p became dramatic with the ad- 
vent of high- mobility Si devices || . Despite the earlier 
attempts to attribute the "anomalous metallic behavior" 
to the density- and temperature-dependent screening || , 
the low-temperature transport in Si MOSFETs fj, as 
well as other 2D systems with a low carrier density n, 
resisted explanation for a decade (for reviews, see J§||). 

Recently, important progress has been made in both 
experiment and theory, which allows to solve this long- 
standing problem. Firstly, it has been recognized that 
multiple valleys in Si MOSFETs enhance interaction ef- 
fects |To|JTl|] . Secondly, the theory of interaction correc- 
tions to the conductivity has been extended beyond the 
diffusive regime pd| ] . This development is crucial because 
the most pronounced increase of p with temperature is 
observed in the ballistic regime (see below) . Thirdly, the 
Fermi-liquid interaction parameters for a dilute electron 
system in Si MOSFETs have been found from detailed 
measurements of Shubnikov-de Haas oscillations |]l2| , |l3| . 
All this allows to compare the experimental data with 
the theory |llf| without any adjustable parameters. 

In this Letter, we show that the most prominent fea- 
tures of the "metallic" state in Si MOSFETs, the strong 
temperature and in-plane magnetic field dependences 
p(T) and p(B^) at relatively small resistivities p <C h/e 2 , 
can be accounted for by the theory |ll[] over a wide 
range of carrier densities, temperatures, and magnetic 



fields. Thus, the "metallic" conductivity of high-mobility 
Si MOSFETs, which was considered anomalous for a 
decade, can now be explained by the interaction effects 
in a system with a large pseudo-spin (two valleys + two 
spins), and an interaction-enhanced spin susceptibility. 
Though the theory accounts only for small correc- 
tions, it works surprisingly well for some samples even at 
relatively high temperatures, where Sp/p ~ 1. The the- 
ory, however, does not agree with experiment at very low 
temperatures and in strong magnetic fields Bu 3> T/gpB, 
as well as at low electron densities, where sample-specific 
effects come into play 

The ac (13 Hz) measurements of the dependences p(T) 
and p(-E>|| ) have been performed on six (100) Si-MOS sam- 
ples from different wafers: Sil5 (peak mobility yP eak = 
4.0m 2 /Vs), Si2Ni (3.4m 2 /Vs), Si22 (3.3m 2 /Vs), Si6-14 
(2.4m 2 /Vs), Si43 (1.96m 2 /Vs), and Si46 (0.15m 2 /Vs); 
more detailed description of the samples can be found in 
Refs. @. 

The quantum corrections to the Drude conductivity 
ud = e 2 nT/m* (r is the momentum relaxation time, m* 
is the effective mass of carriers) can be expressed (in units 
of e 2 /nh) as [11 



a(T, B|| )-a D = 5a c + 15<5ct t + 2[{<r(E z ,T) - o-(0, T)] 
+ 2[a(A v , T) - a(0, T)[ + [a(E z + A v , T) - a(0, T)[ 

+ [a(E z - A v , T) - a(0, T)} + 6a loc (T). (1) 

Here Sa c = x[l - |/(x)[ - ^ ln(^) and Sa T = 
fi(F$)x[l - lt(x, F Q )[ - (1 - f 2 {F§)) £ ln(^) are the 
interaction contributions in the singlet and triplet chan- 
nels, respectively |l7]]; Sa\ oc (T) — ^ ln(r/T v (T)) is the 
weak localization contribution. The terms cr(Z, T) — 
<r(0, T) reduce the triplet contribution when the Zee- 
man energy (Z = E z = 2psB\\), the valley splitting 
(Z = A v ), or combination of these factors (Z = Ez±A v ) 
exceed the temperature. The prefactor 15 to SaT re- 
flects enhancement of the triplet contribution due to 
two valleys of the electron spectrum in (100) Si MOS- 
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FETs 1 10 1 . Because of this enhancement, the "negative" 
correction to the conductivity due to the triplet chan- 
nel, dSax/dT < 0, overwhelms the "positive" correc- 
tion due to the singlet channel and weak localization, 
d(5ac + 5a\ oc )/dT > 0. The triplet-channel prefactor 
might be reduced by inter-valley scattering when the rel- 
evant scattering time t v < h/T. Equation (1) describes 
the quantum corrections in both diffusive and ballistic 
regimes; the crossover occurs at T w (1 + Fq) / (2-kt) [ pd| . 

The terms in Eq. (1) are functions of x = Tt /%, Z, and 
Fq] their explicit expressions are given in Ref. |fu]| . The 
theory Jll| considers the Drude resistivity pu = 1/<Xd 
as a parameter and does not describe the density de- 
pendence /0£)(n). We found pu by extrapolation of the 
linear dependences p(T) to T — 0; r was determined 
from pd using the renormalized effective mass m* in) |l3) . 
The Fermi- liquid parameter Fq = Fq |ll|, which con- 
trols renormalization of the <7*-factor [g* = (75/(1 + Fq), 
where gb = 2 for Si] , has been experimentally determined 
from measurements of the Shubnikov-de Haas (SdH) ef- 
fect 
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Throughout the paper, we accept A„ = 0, 
because small values A„ < IK do not affect the theoret- 
ical curves at intermediate temperatures. Thus, one can 
compare the experiment and the theoretical prediction 
with no adjustable parameters; such comparison is the 
main goal of this paper. In this important aspect, our 
approach differs from the attempts to determine Fq for 
p-type GaAs @ and Si-MOSFETs |oj by fitting p(T) 
and p{B\\) with the theory jLl[ . 
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FIG. 1. Resistivity at B\\ = for sample Si22 vs tem- 
perature. The electron densities, from top to bottom are: 
n = 5.7, 6.3, 6.9, 7.5, 8.1, 8.7, 9.3, 10.5, 11.7, 12.9, 14.1, 
16.5, 18.9, 21.3, 23.7, 28.5, 35.7 (in units of 10 n cm~ 2 ). 
Dots and dashed lines represent the data, solid lines - the 
theoretical curves with A v = and Fg from Ref. Q. 

Figure 1 shows the resistivity of sample Si22 versus 



temperature for different electron densities. A linear 
dependence p(T) extends over a decade in T, up to 
T m O.lEp. In this linear regime, the data almost 
coincide with the theoretical curves (solid lines) calcu- 
lated with no adjustable parameters. Similar agreement 
has been observed for samples Sil5 and Si43. For the 
latter sample, the linear p(T) dependences remain in 
agreement with the theory up to such high temperatures 
(T ~ 0.3E F ) that Sp/p ~ 1 (see Fig. 2). In this case, 
which is beyond the applicability of the theoretical re- 
sults |ll[], we still calculated the corrections to the re- 
sistivity according to Sp — —5op 2 D H). For much more 
disordered sample Si46, the agreement with the theory 
is less impressive: the theoretical p(T) curves are con- 
sistent with the data only at temperatures below 10K, 
which correspond for this sample mostly to the diffusive 
regime. 
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FIG. 2. p(T)-dependences at By = for sample S143 on 
the linear T-scale (main panel) and logarithmic T-scale 
(insert). Dots show the data, solid lines correspond to 
Eq. (1) with A„ = and Fq from Ref. Jl3|. The densities 
are (in units of 10 n cm- 2 ): 1 - 1.49, 2 - 1.67, 3 - 1.85, 4 
- 2.07, 5 - 2.30, 6 - 2.75, 7- 3.19, 8 - 3.64, 9 - 4.54, 10- 
5.43, 11 - 6.33, 12 - 8.13, 13 - 9.91. The border between 
diffusive and ballistic regimes is shown in the inset as the 
dashed line jnj. 

For all samples, the data agree with the theory over a 
broad interval of temperatures (see Figs. 1 and 2), and 
depart from the theory on both sides of this interval. At 
high temperature, the p(T) data for Sil5 and Si22 de- 
viate from the theory up (Fig. 1), while for Si43, they 
deviate mostly down (Fig. 2). The non- universal devia- 
tions indicate the importance of sample-specific localized 
states at T ~ Ep. 

On the low-temperature side, the dependences p(T) 
tend to saturate for all samples, in contrast to the theo- 
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retical prediction Eq. (1). Weakening of the p(T) depen- 
dence might be caused by a non-zero valley splitting at 
temperatures T < A v . Indeed, for samples Si22, Sil5, 
and Si6-14, the saturation temperature (0.2 - 0.5 K, see 
Fig. 1) is of the order of valley splitting estimated from 
SdH measurements, ~ (0.6 — 0.8)K. However, for sam- 
ple Si43, the saturation temperature is too high (1-8 K, 
depending on the density - see Fig. 2), which makes this 
interpretation of the saturation dubious. It is also un- 
likely that the saturation at such high temperatures could 
be caused by electron overheating. One of the reasons 
for decrease in the interaction contribution might be a 
strong (and sample-specific) inter-valley scattering. A 
theory which takes inter-valley scattering into account is 
currently unavailable. 

Finally, we note that for lower densities and higher 
resistivities p ~ h/e 2 , the slope dp/dT changes sign. 
This phenomenon, known as "metal-insulator transition 
in 2D" (2D MIT), can not be accounted for by the theory 
jnj; the corresponding data for samples Si43 and Sil5 
across the 2D MIT can be found in Refs. P,pj[. Here we 
do not discuss this non-universal and sample-dependent 
regime [M. 
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FIG. 3. Magnetoresistivity for sample Si6-14 versus Bl 
at low fields Ez/T < 1 (a), and versus Bu at high fields 
Ez/T > 1 (b). The electron density n = 4.94 x 10 n cm" 2 
is the same for both panels. Lines are the best fits with 
the Fq values shown in the panels. 

We now turn to the magnetoresistance (MR) data. In 
contrast to the temperature dependences of /?, the mag- 
netoresistance agrees with the theory jll|] only qualita- 
tively. For this reason, in fitting the experimental data, 
we treated Ftf as an adjustable parameter. 

Figure 3 shows that the Ffi values, found for sam- 
ple Si6-14 from fitting at low {E Z /2T < 1) and high 
(Ez/2T ;» 1) fields, agree with each other within 10%; 
at the same time, these values differ by 30% from the val- 
ues determined in SdH studies JL3]. A similar situation 



is observed for sample Si43 in weak and moderate fields 
(Fig. 4c, f and Fig. 4b, e, respectively). The weak system- 
atic decrease of |Fq| with B\\ in this field range (Figs. 4) 
agrees qualitatively with non-linearity of magnetization 
which we observed in SdH measurements 

Fitting the weak-field MR at different temperatures 
provides a T-dependent Fq (Fig. 3a). On the other hand, 
our SdH data for the same sample do not confirm such a 
dependence: g*m* is constant within 2-3% over the same 
temperature range |L3| . This discrepancy stems from the 
fact, that, according to our data (see also Ref. B), the ex- 
perimental low-field dependence p(Bu,T) oc B /T a with 
awl, differs from the theoretical one [Eq. (1)]. 
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FIG. 4. Magnetoresistivity for sample Si43 vs Sm and 
B\ at T = 0.3K for two densities: n = 1.49 x 10 n cm" 2 
[panels a, b, c)] and 4.54 x 10 11 cm~ 2 [panels d, e, f)]. The 
upper horizontal scales show P = g* [ibBu/2Ef on panels 
a and d, and Ez/2T on panels b, c, e, and /. 

The discrepancy between the theory and the MR data 
is much more pronounced for sample Si43 in strong fields 
Ez 1 2T S> 1 , where even the sign of deviations becomes 
density-dependent (compare Figs. 4a and 4e). The non- 
universal behavior of the MR has been reported earlier 
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for different samples fijj], and even for the same sam- 
ple cooled down to 4 K at different fixed values of the 
gate voltage |0. The reason for this might be the in- 
teraction of mobile electrons with field-dependent and 
sample-specific localized electron states. 

The Fq(ti) values obtained from fitting the low-field 
MR for three samples are summarized in Fig. 5a. The 
non-monotonic density dependence of Fq and scatter- 
ing of data for different samples indicate that the MR 
is more susceptible to the sample-specific effects than 
p(T). For comparison, we present in Fig. 5 b the Fq- 
values obtained from fitting of the p(T) data for three 
samples, where we treated Fq as a single adjustable pa- 
rameter. In contrast to Fig. 5a, there is an excellent 
agreement between the Fq values extracted from SdH 
measurements and from fitting the p(T) dependences; 
the agreement is observed over a wide density range 
n = (1.5 - 40) x 10 n cm- 2 . 
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FIG. 5. Comparison of Fq (n) values determined from: 
(a) fitting p(Bn) for three samples, and (b) from fitting 
p(T) for three samples. Dashed lines depict upper and 
lower limits for Fq from SdH measurements Jl3|. 

In summary, we performed a quantitative comparison 
of the p(T, n) and p(Bu , n) data with the theory which 
accounts for electron-electron interactions JllJ . For high- 
mobility samples, we found an excellent agreement (with 
no adjustable parameters) between p(T) and the theory 
in the ballistic regime over a wide range of tempera- 
tures and electron densities n = (1.5 — 40) x 10 11 cm~ 2 . 
Our experiments strongly support the theory attribut- 
ing the anomalous "metallic" behavior of high-mobility 
Si MOSFETs to the interaction effects. Non-universal 
(sample-dependent) deviations from the theory pi| have 

5 > 0.1h7P, and 



(b) at low resistivities p <C h/e 2 for both the lowest 
temperatures and high temperatures (T ~ Ep). The 
sample-dependent deviations from the theory are more 
pronounced in the in-plane magnetoresistance, especially 
in high fields (2/Xb-Bm/T > 1). We attribute this non- 
universality to interaction of the mobile electrons with 
field-affected localized electron states. 
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